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ABSTRACT 

We present the age distributions for star clusters and individual stars in the Small Magellanic Cloud 
(SMC) based on data from the Magellanic Clouds Photometric Survey by Zaritsky and collaborators. 
The age distribution of the SMC clusters shows a steep decline, dN duster / dr oc T-O-SSio.is^ ^^gj. ^.j^g 
period 10^ ^ ^ 10^ yr- This decline is essentially identical to that observed previously for more 
massive clusters in the merging Antennae galaxies, and also for lower-mass embedded clusters in the 
solar neighborhood. The SMC cluster age distribution therefore provides additional evidence for the 
rapid disruption of star clusters ( "infant mortality" ) . These disrupted clusters deliver their stars to 
the general field population, implying that the field star age distribution, dNfMstar/dr, should have an 
inverse relation to dNduster/dr if most stars form initially in clusters. We make specific predictions for 
dNfidstar/dr based on our cluster disruption models, and compare them with current data available 
for stars in the SMC. While these data do not extend to sufficiently young ages for a definitive test, 
they are consistent with a scenario wherein most SMC stars formed in clusters. Future analyses of 
dN fidstar / dr that extend down to ages of ~ few million years are needed to verify the age relationship 
between stars residing in clusters and in the field. 
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1. INTRODUCTION 

The age distribution of a population of star clus- 
ters contains important information about the forma- 
tion and disruption of the clusters^. The number of em- 
bedded clusters in the solar neighborhood (with masses 
~ 10^ — 10^ Mq) declines rapidly with age, and it has 
been estimated that only 4 — 7% of these clusters will 
survive for 100 Myr (see the review by Lada & Lada 
2003). We recently found that much more massive (com- 
pact) clusters (> 3x 10^ Mq) in the merging "Antennae" 
galaxies also show a rapid decline in their age distribu- 
tion, with dNciuster/dr oc r"^ (Fall, Chandar, & Whit- 
more 2005, hereafter FCW05). We interpret this decline 
in numbers as due to a high rate of early disruption ( "in- 
fant mortality") for massive clusters in the Antennae. 
Taken together, these age distributions suggest that com- 
pact, massive clusters in the chaotic environment of two 
merging disk galaxies disrupt at approximately the same 
rate as lower mass, embedded clusters in the more quies- 
cent solar neighborhood. This suggests that star clusters 
disrupt rapidly due to processes internal to the clusters 
themselves (e.g., the removal of interstellar material due 
to the mass and energy input from massive stars; Hills 
1980; Fall 2004; FCW05), and implies that the age dis- 
tribution should be roughly dN duster /dr cx t^^ in all 
star-forming galaxies. 

In contrast to the results for clusters in the Milky 
Way and in the Antennae, the cluster age distribution 
in the Small Magellanic Cloud (SMC) has been reported 

^ We use the term "cluster" generically to mean any aggregate 
of stars — regardless of mass, size, or age, and whether bound or 
unbound - with a density significantly higher than the local stellar 
background, making it recognizable as a distinct entity. 



as dNduster/dr oc r^^-^ (Rafelski & Zaritsky 2005; here- 
after RZ05). However, this age distribution is actually 
the number of clusters divided by the number of in- 
dividual stars of the same age, i.e., it is a normalized 
cluster age distribution (RZ05). Does the proposed 
power law for the normalized cluster age distribution in 
the SMC differ from the form found in the Anten- 
nae and solar neighborhood because of the normaliza- 
tion procedure, because clusters in the SMC are forming 
and/or dissolving at a different rate, or for some other 
reason? The motivation for normalizing the number of 
clusters by the number of individual stars was to assess 
whether the formation of stars in clusters parallels the 
formation of stars in the field. An underlying assump- 
tion made by RZ05, then, is that star formation occurs 
in both clusters and in the field simultaneously. How- 
ever, our study of the Antennae galaxies suggests that at 
least 20% and possibly all of the Ha fiux is emitted by 
compact clusters, which disrupt rapidly (FCW05). This 
picture implies instead that at least 20% and possibly all 
stars form in clusters, with stars being subsequently de- 
livered from clusters to the field through the process of 
cluster disruption. This in turn implies that dNduster/dr 
and dN fidstar /dr should have an inverse relationship to 
each other. 

If the SMC cluster system is found to have an un- 
normalized age distribution oc r^^'^, then our picture 
of clusters disrupting at essentially the same rate in all 
star-forming environments may have to be significantly 
revised. What is the unnormalized SMC cluster age dis- 
tribution, does it differ from the normalized cluster age 
distribution, and what is the relationship between field 
and cluster stars? Addressing these questions is the main 
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focus of this Letter. 

2. AGE DISTRIBUTION OF CLUSTERS IN THE SMC 

Wc use the SMC cluster properties published in RZ05 
to directly determine the cluster age distribution. The 
original cluster sample was selected by Hill & Zaritsky 
(2006), and wc refer the reader to their paper for de- 
tails. Briefly, they identified cluster candidates using 
UBVI data from the Magellanic Clouds Photometric 
Survey (MCPS; Zaritsky, Harris, & Thompson 1997). All 
candidate clusters detected in median-subtracted stellar 
images were visually inspected, and only obvious clus- 
ters were retained. Clusters embedded in strong nebular 
emission regions were excluded from the sample due to 
the difficulty of measuring their true integrated colors. 
This leads to a significant bias against clusters younger 
than ^ 10^ yr. Due to these selection criteria and the 
small distance to the SMC, the cluster sample of Hill 
& Zaritsky (2006) is likely to be approximately surface- 
brightness limited, rather than luminosity-limited (as is 
typical for more distant cluster systems). The adopted 
procedure resulted in a catalog of 204 unambiguous SMC 
clusters. 

Basic cluster properties were estimated as follows. 
RZ05 measured integrated U BVI magnitudes within the 
radius enclosing 90% of the V-band light from each clus- 
ter (as tabulated in Hill & Zaritsky 2006). They then 
compared the colors (typically U ~ B, B — V, and V ~ I) 
with the simple stellar population (SSP) model predic- 
tions of STARBURST99 (Leitherer et al. 1999) and 
GALEV (Anders & Pritze-v. Alvensleben 2003) and 
adopted the age corresponding to the minimum x^- RZ05 
provide electronic tables of derived cluster ages and total 
V band huiiinosities for the clusters in their catalog. Wc 
estimated the mass of each cluster by combining the to- 
tal V band luminosity provided by RZ05 with theoretical, 
age-dependent SSP mass-to-light ratios (M/Lv), assum- 
ing a distance modulus of 18.87 (Harries et al. 2003) and 
a foreground reddening value of Eb-v = 0.09. These 
M/Ly ratios assume a Salpeter initial stellar mass func- 
tion (IMF) down to 0.1 Mq, but would be ~ 1.5 times 
lower if we had adopted a Chabrier-type (2003) IMF that 
flattens below 1 Mq . Our procedure does not account for 
any additional extinction within the SMC, which is likely 
higher for younger clusters, resulting in an underestimate 
for the masses of these objects. However, the extinction 
in the SMC is quite small along most lines of sight (e.g., 
Zaritsky 1999), and adopting a single extinction value 
has little effect on our results. 

From the above data we constructed the age distribu- 
tion of star clusters in the SMC, dN duster / dT , by count- 
ing the number of clusters (iV) in different age bins. In 
this Letter, we use cluster ages taken from the RZ05 
comparison with the STARBURST99 model (Leitherer 
et al. 1999) of metaUicity Z=0.004, which is well-suited 
for the young cluster population in the SMC. However, 
we have checked that our basic results are not sensitive to 
the adopted metallicity, population synthesis model, or 
bin size. The la uncertainties are calculated as ±\fN for 
each bin. The data were restricted to cover the age range 
over which the sample of clusters is reasonably complete 
down to ~ 10^ Mq. 

Figure 1 shows the age distributions resulting from 
the full RZ05 sample with log(r/yr) ^ 9.5, and from 
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Fig. 1. Age distributions of star clusters in the SMC and 
the Antennae. The different symbols refer to the approximately 
surface-brightness limited SMC cluster sample of RZ05 (filled cir- 
cles), a mass-limited subset of the RZ05 cluster sample (open cir- 
cles) , and massive star clusters in the Antennae (filled triangles) as 
presented in FCW05. The data point for the youngest bin in the 
SMC cluster age distributions is a lower limit, due to the cluster 
selection procedure (see text). The solid lines represent the power 
law dN^^^ste^r/dT oc r"!. 



the approximately mass-limited sample (> 10^ Mq). 
As expected, the age distribution constructed from the 
mass-limited sample is somewhat shallower than the one 
constructed from the entire cluster sample, since there 
are relatively few clusters more massive than 10^ Mq 
at very young ages. Both samples have a median age 
~ 10^ yr. This is a generous upper limit on the true 
median age due to the original selection, which delib- 
erately excluded clusters in emission-line regions, with 
ages T ^ 10^ yr. Regardless of this incompleteness at 
very young ages, the age distributions of Figure 1 show 
that in the SMC the number of clusters declines mono- 
tonically, with no obvious bends or other features. We 
find that the (unnormalized) cluster age distribution in 
the SMC is well approximated by a power law of the form 
dNciuster/dr cx T--o.85±o.i5^ ^^g^, raugc of ages stud- 
ied ^ T ^ 10^ yr]. This form of the SMC cluster 
age distribution is essentially the same as that found for 
star clusters in the Antennae galaxies (FCW05). Our re- 
sult therefore, provides additional evidence for the rapid 
disruption of clusters due to internal processes ("infant 
mortality"). 

3. FIELD STAR AGE DISTRIBUTION AND RELATIONSHIP 

TO THE CLUSTERS 

This rapid decline of dNduster/dr with r naturally pre- 
dicts that dNfidstar/dr should increase with r as dissolv- 
ing clusters contribute their stars to the field. However 
the exact form of dN fidstar I dr will depend on the frac- 
tion of stars that originally form in clusters. In order 
to explore the relationship between stars in clusters and 
in the field, we use the statistical cluster models pre- 
sented in Whitmore, Chandar, & Fall (2006; hereafter 
WCF06) to predict the field star age distribution. Here, 
we provide a brief summary of the model inputs and 
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assumptions, and refer the reader to WCF06 for details. 
The models allow different cluster formation histories, in- 
cluding a constant rate of formation and gaussian bursts. 
For simplicity, we assume a constant rate of cluster for- 
mation over the life of the SMC, and an initial mass 
function of power-law form with dN duster /dM oc 
as suggested by a number of recent works (e.g., Zhang 

6 Fall 1999; Hunter et al. 2003; WCF06). After clus- 
ters form, they are disrupted in two phases - a rapid, 
mass-independent phase (which is modelled as a power 
law dN/dr oc t^, with index 7), and a longer-term, 
mass-dependent phase that mimics the effects of two- 
body evaporation (modelled as constant mass loss, with 
a rate fi^v = 1 x lO"-'^ Mq yr^^; see FaU & Zhang 2001 
and WCF06 for details). The index of the SMC cluster 
age distribution was found in the previous section to be 

7 = —0.85 ±0.15. Each simulated cluster is composed of 
a number of stars, which are assumed to be drawn from 
a Salpeter IMF with lower and upper stellar mass limits 
of 0.1 Mq and 100 Mq respectively. The only free pa- 
rameter in our model is the fraction of stars, /, that are 
born initially in clusters^. The models track the num- 
bers and ages of stars in the field (i.e., stars which form 
outside of clusters initially or end up there after a cluster 
is disrupted), thus allowing us to predict dNfidstar/dr. 

Figure 2 shows our model predictions of dNfidstar/dr 
for three different values of /: 0.0, 0.5, and 1.0. The 
(solid) curve for / = 1.0 (i.e., all star formation oc- 
curs initially in clusters) shows a steep rise in the num- 
ber of field stars (at very young ages), which is ap- 
proximately the inverse of dNduster/dr. However, by 
^10^ yrs the field star age distribution is essentially flat 
(i.e., dN fidstar I dT cx r''). This effect can be understood 
as follows: when dNdusLer / dr is plotted logarithmically, 
the rate at which the cluster population diminishes is ob- 
served as the slope, 7. However, if dNduster/dr oc r~^'^^ 
were plotted linearly rather than logarithmically, one 
would see that the number of clusters declines quickly 
at first, and then begins to flatten out. The field star 
age distribution is roughly the inverse of the cluster age 
distribution, and so must necessarily rise very quickly as 
stars are dispersed into the field. This increasing func- 
tion essentially "saturates" by ^ 10^ yrs, predicting that 
dNfUsiar/dr varies rapidly only for a short time. 

Is the field star age distribution in the SMC consistent 
with most stars forming in clusters? Harris & Zaritsky 
(2004; hereafter HZ04) present the star formation history 
of the SMC. They use all detected stars (i.e., in both clus- 
ters and in the field), from the MCPS ground-based data. 
However, there is a natural bias in their stellar catalog 
against stars in clusters, due to their inability to pho- 
tometer in dense regions. Therefore, we refer to stars in 
their sample as belonging to the "field". We determine 
the number of SMC stars as a function of age by dividing 
the star formation rates presented in Table 2 of HZ04 by 
the mean stellar mass. Figure 2 shows a comparison of 
the field star age distribution in the SMC with our model 
predictions. This shows that over the observed range the 
field star age distribution is essentially flat (although see 

Our idealized calculations divide stars artificially into two dis- 
tinct categories, i.e., "cluster" and "field". In reality, stars likely 
form in a contiimuni of environments from tightly to more loosely 
clustered aggregates. 




8 

log (T/yr) 

Fig. 2. — Comparison between the SMC field star age distri- 
bution (stars) from HZ04 and predictions for dN fi^^tar/d'^ from 
the WCF06 models (see text for details). The power-law index 7, 
which (luaritifies the rate of cluster disruption, is fixed at —0.85, 
as measured from the SMC cluster age distribution. The only free 
parameter in the model is the fraction of stars, /, that form within 
clusters. Note that even if 100% of stars form initially in clusters, 
the field star distribution flattens by ~ 10^ yrs. 



HZ04 for a detailed analysis of small-scale variations in 
the star formation history) . Unfortunately, the youngest 
age in the HZ04 compilation is just where the predicted 
rise in the field star age distribution begins to flatten, 
making it impossible to determine / with this data set. 
Additionally, the HZ04 data do not discriminate between 
stars in clusters and stars in the field, which results in 
"pollution" of dN fidstar /dr from cluster stars. Never- 
theless, Figure 2 shows that despite dNfi^star/dr being 
rather flat, the SMC field star data do not contradict our 
scenario where most star formation occurs in short-lived 
clusters. Our simulations establish that even if 100% of 
stars formed originally in clusters, their rapid rate of dis- 
ruption creates the field star population so quickly that 
after ^ 10'' yrs the two populations will always appear 
to evolve independently. 

4. DISCUSSION AND CONCLUSIONS 

We have shown here that the (unnormalized) cluster 
age distribution in the SMC has the form dNduster/dr cx 
^-o.85±o.i5 ^j^j that overall the field star age distribution 
is nearly flat (i.e., dNjidsiar/dr oc t°, over the age range 
10'' ^ T ^ lO'' yr). The normalized SMC cluster age 
distribution (the number of clusters divided by the num- 
ber of individual stars of the same age) therefore must 
have the same form as the unnormalized distribution [i.e., 
[dNduster/dr) norm OC r~°-^^='=°-i5] . Our conclusion con- 
tradicts that of RZ05, who found that the normalized 
cluster age distribution is (dNduster/dr) norm oc r~'^'^ 
from the same data. The difference in the power-law 
exponent between our work and that of RZ05 comes 
about because we divided the cluster age distribu- 
tion (dNduster/dr) directly by that for the field stars 
(dN fidstar /dr), whereas RZ05 divided dN duster /dlogr 
by dNfidsiar/dlogr, and then appear to have introduced 
an extraneous factor of t~^. 

Our result that the age distribution of SMC clusters 
has a power-law form dNduster/dr oc t--0-85±o.15 
important implications regarding the physical processes 
that disrupt star clusters. The fact that clusters in vastly 
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different environments (including the chaotic merging 
of two galaxies, a quiescent spiral disk, and a gas-rich 
dwarf galaxy) have age distributions with nearly identical 
forms provides strong evidence that star clusters disrupt 
rapidly due to processes internal to the clusters them- 
selves. Over the first few million years photoionization, 
stellar winds, jets and supcrnovae can remove much of 
the interstellar matter from clusters leaving the member 
stars freely expanding (e.g.. Hills 1980; Boily & Kroupa 
2003a,b). These expanding clusters will be observable for 
~ few X lO'' yrs after becoming unbound, until their sur- 
face brightness becomes so low that they become indis- 
tinguishable from statistical fluctuations of stars in the 
field (FCW05). In addition to the removal of interstel- 
lar material, mass loss from stars themselves will likely 
contribute to the continued unbinding of clusters (e.g., 
Applegatc 1986; Chernoff & Weinberg 1990; Fukushige 
& Heggie 1995). Finally, two-body evaporation will de- 
stroy any remaining clusters on longer timescales (e.g., 
Vesperini 1997; Baumgardt 1998; FaU & Zhang 2001). 
The early disruption ( "infant mortality" ) appears to op- 
erate relatively independent of cluster mass, while two- 
body evaporation is mass-dependent. 

The stars from disrupted clusters add to the general 
field, requiring a close relationship between the age dis- 
tributions of clusters and field stars if most stars form 
initially in clusters. In this work we compared age distri- 
butions for clusters and individual stars in the SMC with 
model predictions to further test this general framework. 
Our models suggest that the rapid drop in the cluster age 
distribution should be reflected in a rapid rise in the field 
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star age distribution. However this rapid rise in the num- 
ber of field stars is only observable for the first ^ 10^ yrs 
regardless of the fraction of stars that initially formed in 
clusters. Because the present data on SMC stars begin 
just as the predicted field star distributions flatten, and 
because there is no discrimination between stars belong- 
ing to clusters versus to the field, they do not provide an 
unambiguous answer. However the observed relations are 
consistent with a scenario where most stars form in clus- 
ters that are subsequently disrupted by the mechanisms 
discussed above. The models also predict that even if 
all stars form initially in clusters, galaxies will contain 
a relatively small number of observable clusters at any 
given time. Future studies with the Hubble Space Tele- 
scope (which allow better discrimination between stars 
in clusters and in the general field), and in the infrared 
(which can probe very young, deeply embedded stars and 
clusters) could be used to search for the expected rise in 
the field star age distribution at r ^ 10^ yrs. A measure- 
ment of the rise in dNjidsLar/dr would provide a direct 
estimate of the fraction of stars which originally form in 
clusters versus in the field. 
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